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ABSTRACT 

Context. The enigmatic object HD 45166 is a qWR star in a binary system with an orbital period of 1.596 day, and presents a 
rich emission-line spectrum in addition to absorption lines from the companion star (B7 V). As the system inclination is very small 
(i = 0.77° + 0.09°), HD 45166 is an ideal laboratory for wind-structure studies. 

Aims. The goal of the present paper is to determine the fundamental stellar and wind parameters of the qWR star. 
Methods. A radiative transfer model for the wind and photosphere of the qWR star was calculated using the non-LTE code CMFGEN. 
The wind asymmetry was also analyzed using a recently-developed version of CMFGEN to compute the emerging spectrum in two- 
dimensional geometry. The temporal-variance spectrum (TVS) was calculated for studying the line-profile variations. 
Results. Abundances, stellar and wind parameters of the qWR star were obtained. The qWR star has an effective temperature of 
r e jf = 50000 ± 2000 K, a luminosity of log(L/L Q ) = 3.75 ± 0.08, and a corresponding photospheric radius of ^ p hot = 1-00 Rq. 
The star is helium-rich (N(H)/N(He) = 2.0), while the CNO abundances are anomalous when compared either to solar values, to 
planetary nebulae, or to WR stars. The mass-loss rate is M = 2.2 x 10~ 7 M Q yr~' , and the wind terminal velocity is v„ = 425 km s . 
The comparison between the observed line profiles and models computed under different latitude-dependent wind densities strongly 
suggests the presence of an oblate wind density enhancement, with a density contrast of at least 8:1 from equator to pole. If a high 
velocity polar wind is present (~ 1200 km s~'), the minimum density contrast is reduced to 4: 1. 

Conclusions. The wind parameters determined are unusual when compared to O-type stars or to typical WR stars. While for WR 
stars v„/v esc > 1.5, in the case of HD 45166 it is much smaller (v^/v^ = 0.32). In addition, the efficiency of momentum transfer 
is j] = 0.74, which is at least 4 times smaller than in a typical WR. We find evidence for the presence of a wind compression zone, 
since the equatorial wind density is significantly larger when compared to the polar wind. The TVS supports the presence of such a 
latitude-dependent wind and a variable absorption/scattering gas near the equator. 

Key words. Stars: winds, outflows - Stars: mass-loss - Stars: fundamental parameters - binaries: spectroscopic - Stars: individual: 
HD 45 166 - Stars: Wolf-Rayet 



5_i 1 . Introduction spectrum of HD 45166. By fitting the continuum energy distri- 

. .. . bution from the far ultraviolet to the near-infrared, they derived 

HD 45166 has been observed since 1922 (|Anger||l933|), with- that the qWR star has R = \ 2 Q R^ log(L/L ) = 3.28, and 

out much ad vancement in the understanding of its nature. TeS = 35000 K A maS s-loss rate of 2 x 1(T 7 M^r 1 was in- 

|van Blerkoml U978J) analyzed the optical H and He lines as- f erre d by those authors using the optical lines of He n, C m, and 
suming that HD 45166 is a Population I WR object, and con- 

eluded that the WR component has a radius of 1 R , and has a ISteiner & Olivekal (120051) (hereafter Paper I) showed that 

small-sized envelope expanding with a velocity of 150 km s 1 HD 45166 is a doub i e -lined binary system composed of a qWR 

The resulting number density of He II is about 10 u cm- and a B7 v star in a system with an orbital period of j 5% day 

and mimics the environment of a WR envelope. He found a HD 45166 presents a rich emission-line spectrum in addition to 

mass-loss rate of 4.5 x 10 s M yr and wind densities of the absorption spe ctrum due to the cooler component. The or- 

N(He) = 3 x 10" cm- J and WH)=10 12 cm J . Analyzing the bital parameters of the system were derived in Paper I, showing 

IUE data, [Willis & Stickland j 09831) derived log(L/L ) = 3.84, that the orb i t is slightly eccentric (e = 0.18 ± 0.08) and has a 

a radius of R = 0.77 R^_and an effective temperature of very small inc ii na tion angle (/ = 0.77° + 0.09°). The masses are 

r eff =60000 K. [Willis etal. | dl989j) obtained a wind terminal ve- M[ = 4 2 Mq and Ml = 4 8 Mq . In addition to the orbital period, 

locity of 1200 kms ', derivedfrom the UV resonance lines. two other periods were found in me qWR star (5 an d 15 hours, 

Recently, | WiUis & Burnleyl d2006j) analyzed the far-ultraviolet p aper I} As me system i nc ii na tion is very small, HD 45166 is 

an interesting laboratory for studying the wind structure. 

* Based on observations made with the 1.52 m ESO telescope at La The goal of the present paper is to study in detail the stel- 

Silla, Chile. far and wind parameters of the qWR star. For this purpose 

** e-mail: jgroh@mpifr-bonn.mpg.de we use the radiative transfer code CMFGEN (Hillier & Miller 
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1998) to analyze the high-resolution optical spectrum presented 
in Paper I. The tem poral variance spectrum (hereafter TVS; 
iFullerton et al.|[l996l) of the strongest emission lines is also ana- 
lyzed in order to obtain insights on the wind structure. 

This paper is organized as follows. In Sect. [2] we briefly sum- 
marize the data presented in Paper I, and describe how the spec- 
trum of the qWR star was disentangled from the B7 V compan- 
ion. In Sect. [3] we present the main characteristics of CMFGEN, 
while in Sect. [4] the results of the quantitative analysis using 
the spherical modeling are shown. The latitude dependence of 
the wind is analyz ed in Sect. [5] using a rec ently-developed ver- 
sion of CMFGEN dBusche & Hillieril2005l) . The analysis of the 
TVS is presented in Sect. [6] and the results obtained in this 
work are discussed in Sect. [7] especially the presence of a wind- 
compression zone. Finally, Sect.[8]summarizes the main conclu- 
sions of this paper. 

2. Observations and subtraction of the B7 V 
companion spectrum 

The data analyzed on this paper were taken in January 2004 with 
the Fiber-fed Exte nded Range Optical Spectrograph (FEROS, 
Kauf er et al] [T999) at the 1.52 m telescope of the European 
Southern Observatory (ESO) in La Silla, Chile. The spectra have 
a resolution power of R=48 000, and were red uced using the 
standard data-reduction pipeline (IStahl et all 19991) . A total of 40 
spectra with individual exposure times of 1 5 min were averaged, 
resulting in a total integration time of 10 h. 

Since our goal in this paper is to obtain the fundamental pa- 
rameters of the qWR star, its spectrum has first to be disentan- 
gled from the spectrum of the B7 V companion. This task was 
accomplished by: 

1 . flux calibrating the average observed spectrum of HD 45 166 
using the available photometry (de-reddened); 

2. scaling the flux of a s tandard synthe tic B7 V continuum spec- 
trum, obtained from Pickles ( 1998), to a distance of 1.3 kpc 
(Paper I); 

3. multiplying the normalized spectrum of the B7 V star 
HD 62 714 obtained from t he UVES Paranal Observatory 
Project (Bagnu lo et al.l2003h by the Pickles B7 V scaled con- 
tinuum obtained above, to account for the spectral features of 
the B7 V companion; 

4. subtracting this scaled B7 V spectrum from the flux- 
calibrated observed spectrum of HD 45166; 

5. normalizing the resultant spectrum by the continuum, which 
can then be directly compared to the continuum-normalized 
model spectrum. 

In general, this procedure results in a very good subtraction 
of the continuum and reasonable subtraction of the hydrogen ab- 
sorption lines due to the B7 V star, as can be seen in Fig. Q] 
However, the rest of the rich absorption-line spectrum of the 
companion star is not subtracted perfectly, and weak residual 
lines can still be seen. Nevertheless, their presence do not af- 
fect the diagnostic lines used for obtaining the parameters of the 
qWR, as the companion lines are quite narrow, weak, and often 
are not blended with the lines of the qWR star. 

3. The Model 

We u sed the radiative transfer code CMFGEN dHillier & Milled 
Il998l) to analyze the spectrum of the hot component of 
HD 45166 in detail. The code assumes a spherically-symmetric 



3xl0" 12 " 




X 

E i.5xio" r 



l.Oxio" 12 r 




7500 8000 8500 

Wavelength (A) 

Fig. 1. Resulting de-reddened spectrum of HD 45166 decom- 
posed into the qWR star (full black line) and the B7 V compan- 
ion (dashed red line). As anticipated in Paper I, the qWR star 
is dominant in the blue part of the spectrum, while the B7 V 
companion dominates for A > 5000A. 



outflow in steady-state, computing line and continuum forma- 
tion in the non-LTE regime. Each model is specified by the hy- 
drostatic core radius the luminosity L*, the mass-loss rate M, 
the wind terminal velocity v m , and the chemical abundances Z, of 
the included species. Since the code does not solve the momen- 
tum equation of the wind, a velocity law must be adopted. The 
velocity structure v(r) is parameterized by a beta-type law, which 
is modified at depth to smoothly match a hydrostatic structure 
at R+ (defined at a Rosseland optical depth of r = 100 in our 
models). The use of a hydrostatic structure at depth is crucial to 
analyze HD 45166, since the continuum and some of the higher 
ionization lines are formed in this region. 

The code includes the effects of clumping via a volume fill- 
ing factor / which depends on the distance following an expo- 
nential law that supposes an unclumped wind when r — R+. The 
clumps start to form at the distance where v = v c , and the wind 
reaches its maximum clumping at large r: 

/(r) = / + (l-/)exp[-v(r)/v c ]. (1) 

CMFGEN includes directly the influence of line blanket- 
ing in the ionization structure of the wind and, consequently, 
in the calculated spectrum. With the concept of super-levels, the 
equations of statistical equilibrium and radiative transfer can be 
solved simultaneously including thousands of lines in NLTE. 
The atomic model used for HD 45166 consists of 35059 spectral 
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Table 1. Final atomic model used in the analysis of HD 45 166. A 
single-level ion of one higher ionization stage is included for ev- 
ery species (e.g. C vi for carbon) and omitted below for brevity. 



Ion 


N s 


N F 


N T 


Hi 


30 


20 


435 


He i 


45 


40 


435 


He ii 


30 


22 


233 


Cm 


243 


83 


5513 


Civ 


64 


60 


1446 


C v 


107 


43 


2196 


Nm 


287 


57 


344 


Niv 


60 


34 


331 


Nv 


67 


45 


1554 


Om 


349 


64 


823 


Oiv 


95 


19 


823 


O v 


34 


34 


118 


Ovi 


13 


13 


59 


Si iv 


33 


22 


185 


Siv 


45 


22 


302 


Feiv 


264 


28 


3889 


Fe v 


191 


47 


2342 


Fe vi 


433 


44 


8663 


Fe vn 


252 


41 


2734 


Fe vm 


324 


53 


7305 


Niiv 


200 


36 


2337 


Niv 


183 


46 


1524 



TVs = number of included superlevels; Np = number of total ener- 
gy levels; Nt = number of bound-bound transitions included for each 
ion. 



lines from 2426 energy levels grouped in 617 super- levels of H, 
He, C, N, O, Si, Fe, and Ni, as shown in Table Q] 



4. Results: fundamental parameters of the qWR star 

4.1. Effective temperature, luminosity and radius 

The effective temperature was constrained using the relative 
strength of lines from different ionization stages of He, C, and 
N. The diagnostic lines used to obtain the He ionization struc- 
ture were Hei A 5876, Hei A 6678, Hen A 4686, and Hen A 
541 1, while for the C ionization structure we used C in AA 4647- 
4650^1651 and Civ AA 5801-5812. The N ionization structure 
was obtained using the lines of N m A 4097, N m A 4634, N iv 
A 4058, Niv AA 7103-7109-7123, Nv A 4605, Nv A 4620, and 
NvA4945. 

Figure |2]presents the fits to the diagnostic lines used to de- 
rive T e ff in this work. In particular, the relative strength between 
He i and He n lines requires T e g > 47 000 K, otherwise the He i 
lines in the model become too strong compared to the obser- 
vations, and the He n lines become too weak. The N ionization 
structure requires models with T e g > 46 000 K in order to re- 
produce the ratio between the Nm and Niv lines mentioned 
above. Reasonable fits to the N v lines are only achieved by 
models with T^g > 50 000 K. The strongest C lines seen in the 
spectrum, namely C iv AA 5801-5812, also require models with 
Teff > 49 000 K to obtain reasonable fits to the observations. 
CMFGEN models with T eff < 48 000 K produce too strong Cm 
emission. 

As can be seen in Fig. [2] the fits to the observed line spec- 
trum are very sensitive to the effective temperature, and the use 
of three diagnostics in this temperature regime allowed us to 
constrain the effective temperature as r e jf(T = 2/3) = 50000 + 
2000 K. The temperature of the hydrostatic core was constrained 




m J , ; I |- , | | | | , | | | | , ^ | , , . , | 1 H ■ ■ , , I ■ , , ■ I , . 

4941 4946 5790 5805 5820 5870 5880 7100 7110 7120 



Wavelength (A) 

Fig. 2. Normalized spectrum of the qWR star in HD 45166 
(black line) compared with spherical CMFGEN models with 

- 50000 K (purple dash-double-dotted line), T e g 
46 000 K (blue dot-double-dashed line), T ef{ = 43 000 K (green 
dot-dashed line), and T e g = 35 000 K ( red dashed line). A mode l 
with the same parameters proposed by Willis & Burnlevl |2006) 
(Teff = 35 000 K, Voo = 1200 kms _I ) produces fits similar as 
the original T e ff = 35 000 K model, and does not fit our optical 
data as well. For clarity, such a model is compared only to Hen 
A 4686 (brown thick line). 



to J*(t = 100) = 70 000 + 2000 Kand obtained through the use 
of a hydrostatic structure at depth. The hydrostatic structure is 
sensitive to the effective gravity and, hence, to the adopted stel- 
lar mass (M = 4.2 M Q , Paper I). 

It is worthwhile noting that even if it was not possible to 
reproduce lines corresponding to the observed ionization stages 
of He, C, and N using a unique best model with T e g = 50 000 K, 
small changes in the range +2000 K were sufficient to adjust the 
discrepant lines (Fig.|2]i. 

Models with other parameter regimes do not fit the opti- 
cal spectrum of HD 45166 from 2004. For instance, increas- 
ing the terminal velocity of a model with = 35 000 K 
to Voo - 1200 km s" 1 (i.e., the same model parameters of 
Willi s & Burnl ey 2006) does not enhance the ionization degree 
of the wind and, therefore, does not fit the high-ionization opti- 
cal lines. Such high-terminal velocity model provides fits similar 
to the original Teff = 35 000 K model with v m = 425 km s _1 (Fig. 

The fundamental parameters of the qWR star obtained in this 
work, using tailored CMFGEN models to fit the continuum en- 
ergy distribution and the optical spectrum, differ significantly 
from the previous published values. While we cannot exclude 
variability, the discrepancies between our results and previous 
works are at least partially explained by the inclusion of key 
physical ingredients in our analysis; namely, the inclusion of de- 
tailed non-LTE radiative transfer in the co-moving frame, full 
line blanketing, and use of a hydrostatic structure at depth, which 
allows simultaneous modeling of the photosphere and wind. 
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However, since Wil lis & Burnlevl ([2006) used the same radiative 
transfer code as ours, the aforementioned physical ingredients 
were presumably included in their modeling as well. The value 
of T eS inferred by those authors (35 000 K) is 15 000 K lower 
than our determination, and this difference is puzzling. 

In order to investigate this discrepancy in the values of T e ff, 
we computed several CMFGEN models covering a broad range 
of effective temperatures in order to analyze the changes in the 
continuum, sin ce this was the diagnost ic of the effective temper- 
ature used by Willi s & Burnle y (2006). While the spectral lines 
are strongly sensitive to changes in T e ff (Fig. we found that 
the continuum slope predicted by CMFGEN models in the range 
32 000 K < 7; ff < 50 000 K is almost insensitive to the adopted 
value of Teff (Fig. [3^,b). Indeed, all such models provide good 
fits to the ultraviolet-to-near-infrared SED of HD 45166 (Fig. 
[3};). This is not surprising since such hot stars emit the bulk of 
their flux in the range 228-1000 A. Similar results ha ve long 
been found for other hot star s, such as Wolf-Ra vet stars (iHillierl 
19871; lAbbott&Contilll987l) and O-type stars dMartins & Plezl 
2006). We also computed models using a large value of the wind 
terminal velocity, as proposed by Willis & Burnley (2006), and 
found no noticeable changes in the continuum slope. This hap- 
pens because the photosphere is located at rather low velocities 
(~ 30 km s _1 ) in comparison with Wolf-Rayet stars, and changes 
in the velocity field have little impact on the continuum forma- 
tion region and in the continuum spectrum shown in Fig. [3] 

In the case of HD 45 166, only small changes (of the order of 
8% or less) are seen in the ultraviolet when comparing models in 
the range 32000 K < T eff < 50000 K (Fig.E>). Those changes 
are easily compensated by slightly adjusting other model param- 
eters, such as the mass-loss rate, core radius, or reddening law. 
Even without changing any parameter, those small changes in 
the ultraviolet flux tend to be masked by observational errors, 
since the lcr error in the absolute photometry of the IUE and 
FUSE flux-calibrated spectrum is at least 5%. Therefore, since 
models in the range 32 000 K < T e g < 50 000 K provide good 
fits to the observed SED (Fig.\3]p), we suggest that the technique 
of fitting the continuum is not suitable for constraining T e g in 
the case ofHD 45166, and very li kely explains the discrep ancy 
between the low T e $ suggested by \Willis & Burnle^ {2006) and 
the value determined in this work. 

The stellar luminosity of the qWR star was obtained by 
matching the best-model flux, scaled to a distance of d — 1.3 kpc 
(Paper I), with the observed flux of HD 45166, de-r eddened us- 
ing E( B—V)—0. 155. The photometry was taken from Willi s et alj 
d 19891) . from Paper I, and references therein. We obtained a lu- 
minosity of Li, = 5650 L (log(L/L ) = 3.75), with an error 
due to the uncertainties in the modeling, distance, reddening law, 
and photometry amounting to about 20% (0.08 dex). Figure [3}; 
displays the de-reddened observed flux compared with the flux 
from the best CMFGEN model. We determined a luminosit y 
significantly higher than that given bv lWillis & Burnlevl d2006l) . 
which is due to the significantly higher T e s inferred by our de- 
tailed modeling of the spectral lines. 

Combining the derived values of T+ and T e g with L*, and 
using the Stefan-Boltzmann law, it is possible to determine the 
radius of the hydrostatic core of the qWR star as /?*(t = 100) = 
0.51 R , and the radius of the photosphere (defined as r = 2/3) 
asfl phot = 1.00 Rq. 




i i i i i i i i i i i 

1000 2000 3000 5000 10000 20000 



Wavelength (A) 

Fig. 3. a) Comparison between the continuum flux predicted by 
CMFGEN models for the qWR star in the range 32000K < 
r e ff < 50000K. All models were scaled to the observed de- 
reddened flux of the qWR star at 5500 A (3.09 x 106-13 erg 
s _1 cirT 2 A _1 , see Fig. [1). Since the continuum slope is very 
weakly sensitive on the value of Jeff (see text for discussion), 
the models are basically superimposed to each other, b) Ratio 
between the continuum flux predicted by a CMFGEN model 
with a given r e ff and the continuum flux of the T e ff = 50 000 K 
CMFGEN model, c) Comparison between the total de-reddened 
observed flux of HD 45166 (black) with the flux predicted by 
the radiative-transfer model (red dashed line). The model shown 
here is composed by the sum of the continuum of the B7 V tem- 
plate from the companion (shown in Fig. [TJ and the continuum 
of the best CMFGEN model for the qWR star. 

4.2. The mass-loss rate 

We derived M = 2.2 x 10~ 7 M yr~' in order to reproduce 
the strength of the spectral lines present in the optical spec- 
trum taken in 2004. We also obtained a volume filling factor of 
/ = 0.5, using the electron-scattering wings of Hen A 4686 and 
Civ AA 5801-5812 as diagnostics. However, the electron scat- 
tering wings are weak, and we cannot rule out the presence of 
an unclumped wind (/ = 1), with M = 3.1 x 10~ 7 M yr _I . The 
value of the clumped M derived in this w ork agrees with the pro- 
posed value given bv lWillis & Burnle v (2006), while the value 
of the unclumped M derived here is 50% higher than the value 
obtained by them. 

Using the value of the mass-loss rate combined with the ve- 
locity law derived in Sect l4.3l and the equation of mass continu- 
ity, we obtain the wind density structure shown in Fig. [4] 

4.3. Wind terminal velocity and momentum transfer 

From our optical spectrum we did not detect velocities above 
+550 km s -1 for the He n A 4686 line (see Paper I). For the He i 
A 5876 line, however, the maximum velocity observed is much 
smaller (280 km s" 1 ). For the Civ AA 5801-5812 lines, a strong 
emission exists at low velocity, while weak emission can be seen 
extending up to 600 kms . From the model fit to Hen A 4686 
and He n+Ha A 6560 , we obtained a wind terminal velocity of 
= 425 + 50 km s _1 . The value of the acceleration parameter B 
was set to 4.0 in order to reproduce the relative strength among 
the He n lines of the Pickering series. In particular, it was impos- 
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Table 2. Fundamental parameters derived for the qWR star in 
HD 45166. 



Parameter 



Value 



log (U/L Q ) 


3.75 ± 0.08 


T* (K) 


70000 ± 2000 


T eB (K) 


50000 ± 2000 


R,/R Q (t = 100) 


0.51 


tfpho./Ro (t = 0.67) 


1.00 


M(M Q yr') 


2.2 x 10~ 7 


v„ (km s~') 


425 


P 


4.0 


f 


0.5 


v c (kms ') 


100 



sible to reproduce the high-order He n lines of the Pickering se- 
ries using lower values of ft. The velocity law obtained is shown 
in Fig.|4](panel a). 

The value of v^, derived from the optical lines in this work 
is ~ 3 times lower than what was derived fr om the observations 
of UV resonance lines (v m ~ 1200km s" 1 . IWillis & Stic kland 
1 1 983b . One possible explanation for the different values of Voo 
derived from the UV and optical spectrum is the presence of a 
latitude-dependent wind, with a higher wind terminal velocity in 
the polar direction, and a slower equatorial wind. This possibility 
is further explored in Sect. [5] 

From the optical lines, we obtained that the efficiency of 
the momentum transfer from the radiation to the gas (77) is 
much smaller in HD 45166 than in WR stars. The momen- 
tum of the gas is Mv^, while the momentum of the radia- 
tion is L*/c. T he ratio between the two momenta for WN stars 
is r] = 2.8 dCrowtheri [2007b . while for HD 45166 we ob- 
tain j] = 0.74 . As a consequence, the wind driven by radia- 
tion pressure is less efficient in HD 45166 than in WR stars, 
which can explain the low ratio Voo/vesc found for HD 45166. 
For instance, O-type stars in about the same temperature range 
have Vco/Vesc = 2.6 dLamers et al.lll995l). while WR stars have 
1.5 < v ra /v esc < 4 (Lamers & Cassinelli 1999). However, in 
HD 45166, v esc =1320 kms -1 , v M =425 kms -1 , and therefore 
Voo/vesc = 0.32, i.e, at least 5 times smaller than in typical WRs. 
This is difficult to achieve for a normal wind, since it requires a 
fine balance between gravity and line forces. The presence of the 
close B7 V companion, fast rotation, and/or deviations from a 
spherical wind could explain such an extrem ely low ratio. Using 
Voo=1200 kms^derived from the UV fines dWillis & Stick land 



1983), the ratio v^/v^ is increased to ~ 1, becoming closer to 
the value for O-type and WR stars. 




log (r/R,) 



log (r/R.) 



Fig. 4. Velocity law (panel a), density structure (b), electron 
density (c), and Rosseland optical depth (d) of HD 45166 as a 
function of distance. 



Table 3. Surface chemical abundances of HD 45166. 



Species Number fraction Mass fraction Z/Z e 



H 


2.0 


3.3x10-' 


0.46 


He 


1.0 


6.5x10' 


2.34 


C 


3.0X10- 3 


5.9xl0- 3 


1.93 


N 


2.0xl0~ 3 


4.6xl0~ 3 


4.17 


O 


1.5xl0- 3 


3.9xl0 3 


0.41 


Si 


6.8xl0 -5 


3.6xl0 -4 


1.00 


Fc 


1.3xl0 -4 


1.4xl0 -3 


1.00 


Ni 


6.6xl0 -6 


7.3xl0 -5 


1.00 



Table 4. Abundance ratios (in number) of HD 45166 compared 
to different object classes. 



Object 


He/H 


C/N 


O/N 


Reference 


HD 45166 


0.5 


1.5 


0.75 


1 


Sun 


0.095 


8.6 


7.2 


2 


PNI 


0.14 


0.7-3 


1.2 


3 


PNII 


0.11 


1.4-9 


3.5 


3 


WNL 


H-free 


0.04 




4 


WNE 


H-free 


0.02 




4 


WNEw 


1.5 


0.005 


0.05 


4 



( 1990); (4) ICrowthej (2007). 



4.4. Surface abundances 

The surface abundances obtained with CMFGEN are summa- 
rized in Table [3] The helium abundance was determined adjust- 
ing the helium content in the models to match the observed in- 
tensity of He n A 4686, He 11 A 541 1 . and the blend of He ii+Hq- 
A 6560. The value obtained is N(H)/N(He)=2.0, which corre- 
sponds to about 2.3 times the solar abunda nce of He (hereafter 
we used the solar abundance values from ICoxl 12000 and refer- 
ences therein). 

The model calculations constrained C and N abundances that 
are significantly higher than the abundances of O and Si. The 
abundances of HD 45 166 are in fact quite anomalous when com- 
pared to solar values, or to the values found in central stars of 
planetary nebulae or Wolf-Rayet stars (see Table [4|. This em- 
phasizes the peculiar nature of HD 45166. 



The chemical abundances must be analyzed in view of the 
evolutionary status of both stars in the HD 45166 system, which 
is beyond the scope of this paper. Such a detailed analysis will be 
the subject of a forthcoming Paper III. We antici pate that the star 
is likely an exposed He core (Marco et al. 2007), which is prob- 
ably related to the presence of the close secondary. We suggest 
that He burning is currently going on in the nucleus, producing 
carbon, and, thus, explaining the He and C overabundance. As 
the star still has an H-rich envelope, there certainly is a shell 
burning H through the CNO-cycle, enhancing the N content. As 
can be seen in Table [3] heavier elements such as Si, Ni and Fe 
have solar abundances. 
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Fig. 5. Ionization fractions of HD 45166 as a function of dis- 
tance, derived from the spherical CMFGEN models of the most 
abundant ions in the wind of the qWR star. The ionization frac- 
tions are normalized by the total number of ions. 



4.5. The ionization structure of the wind 

Figure[5]presents the ionization structure of the wind of the qWR 
star determined using the spherical CMFGEN model. H is fully 
ionized along the whole wind, but the He, C, and Fe ionization 
structures are more stratified, with He 2+ and C 4+ being more 
abundant in the photosphere and in the inner wind, while He + 
and C 3+ dominating at distances larger than about 10 R*. This 
also explains the high sensitiveness of the He, C, and Fe lines on 
the model parameters. On the other hand, the O and N ioniza- 
tion structures are dominated respectively by 3+ and N 3+ ions 
in most part of the wind - 4+ and N 4+ are present in significant 
fraction only at distances smaller than ~ 1.5 R*. 
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Fig. 6. Comparison between the line profiles predicted by 
the best CMFGEN spherical model with the observations of 
HD 45166, in the spectral region 3850-5950 A. The strongest 
spectral lines used in the analysis are identified. 



5. Presence of a latitude-dependent wind 

5.1. Observational evidence from the line profiles 

The spherical models obtained with CMFGEN reproduce the 
strength of most of the emission lines, provide a superb fit to 
the continuum of the qWR, and a reasonable fit to the spec- 
tral lines, considering the complex physical nature of HD 45166 
(see Figs. [6] and [7}. However, some discrepancies are present 
when comparing the line profiles predicted by the best spherical 
model with the observations (Fig. [2j». These discrepancies might 
provide key insights on the validity of the model assumptions. 

As already mentioned in Paper I, the hydrogen and he- 
lium line profiles are clearly different from the CNO lines pro- 
files. The later can be very well fitted by Lorentzian profiles, 
while Hen A 4686 and the other Hen lines are better fitted 
by a Voigt/Gaussian profile. The full widths at half maximum 
(FWHM) are also significantly different between the two groups 
of lines. 

Insights on the physical interpretation of the line profiles 
from different species can be obtained with the best CMFGEN 
model obtained in Sect.|H which assumes a spherical wind with 
a monotonic velocity law. We present in Fig. [8] the line for- 
mation regions for the most important diagnostic lines present 
in the spectrum of HD 45166. The inner layers correspond to 
the formation region of higher-ionization lines such as N v and 
C iv, while the outer layers correspond to the formation region of 
lower-ionization lines such as N m, C m, and He i. As can be seen 



in the CMFGEN model spectrum displayed in Figs. [2] [6] and|7] 
the higher-ionization lines are predicted to be narrower than the 
lines of lower ionization. In addition, recombination lines such 
as those from H and He are emitted in an extended region of the 
wind and, therefore, should be broader than N iv and N v lines, 
for instance. 

Indeed, the N iv and N v lines are narrow in the observations, 
while the He n lines are broad, and both are well reproduced by 
the spherical CMFGEN models. However, the lower-ionization 
lines, such as C m, N m, and He i, are surprisingly narrow in the 
observations - much narrower than the He n and C iv lines. 

We examined the effects of changing one or more physical 
parameters to improve the fit of the lower ionization lines, run- 
ning a large number of additional CMFGEN models. None of 
them improved the fit to those line profiles, i.e., they did not 
provide narrower lower-ionization lines. Actually, the only way 
to change the width of those lines using a spherical model is 
using a much lower terminal velocity of ~ 70 km s _1 , which ob- 
viously did not fit the other strong spectral lines. Moreover, it 
would be hard to explain the origin of such a very low wind ter- 
minal velocity in a very hot star (7efr=50000 K). In summary, 
a spherical symmetric wind cannot reproduce the profile of the 
lower-ionization lines. 

Therefore, we propose the presence of a latitude-dependent 
wind in the qWR star to explain the narrow profiles observed in 
the N m, C in, and He i lines. Using the code outlined in Sect. 15. 21 
we examine the effects due to changes in the wind density and 
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spherical stellar wind as in HD 45166. However, this is cur- 
rently impossible to be done taking into account all the relevant 
physical processes, full line blanketing, and the degree of de- 
tail achieved by spherically-symmetric codes such as CMFGEN. 
The large parameter space to be explored in 2D models, and the 
huge amount of computational effort demanded, do not allow to 
perform a self-consistent 2D analysis of very complex objects 
such as HD 45166. 

Nevertheless, it is important to note that significant pro- 
gresses ha ve been achieved in developing full 2D rad iative trans- 
fer codes dZsargo et al .1 l2006t iGeorgievet all 12006). which can 
be applied to HD 45166 in future works. Therefore, it is desir- 
able to explore the parameter space suitable for HD 45 166, using 
justified assumptions to make the computation of the observed 
spectrum in 2D geometry a tractable problem. 

In this work we used a recentl y-developed modification in 
CMFGEN dBusche & Hillierll2005b to compute the spectrum in 
2D geometry. We refer the reader to that paper for further details 
about the co de, whose main characte ristics are outlined below. 

With the iBusche & Hillierl ([2005) code, it is possible to ex- 
amine the effects of a density enhancement and changes in 
the velocity field of the wind as a function of latitude. The 
code uses as input the ionization structure, energy-level popula- 
tions, temperature structure, and radiation field in the co-moving 
frame, as calculated by the original, spherically-symmetric 
CMFGEN model. Usi ng CMF_FLUX (lHilher & Milled 1 1998b 
Busc he & Hilli er 2005), the emissivities, opacities, and specific 
intensity J are calculated from the spherically-symmetric model. 

A density enhancement and wind terminal velocity variation 
can then be implemented using an arbitrary latitude-dependent 
density/wind terminal velocity distribution. We examined the ef- 
fects due to oblate and prolate density parameterizations, which 
have respectively the form 



P2D * p lD (l + a.cos 6) 



(2) 



where 9 is the latitude angle (0°=pole, 90°=equator). The 
changes in the wind terminal velocity as a function of lati- 
tude were parameterized using three parameters ( VB 1 , VB2, and 
VB3), 



v(r) 2D /v(r) 1D = VBl + (VB2 



\cos(6) VB3 \) 



(3) 



We assumed in this work a scaling law according to which 
the 2D model has the same mass-loss rate as the spherically- 
symmetric ID model. 

After computing the density changes and velocity field vari- 
ations, the 2D source function, emissivity, and opacity are cal- 
culated, assuming that these quantities depend only on the new 
values of the scaled density. Consistent scaling laws are used 
for different processes (e.g. density-squared scaling for free-free 
and bound-free transitions, and linear-density scaling for elec- 
tron scattering). 

In the final step, the code computes the spectrum in the ob- 
server's frame, which can then be compared with the observa- 
tions. 



terminal velocity as a function of the stellar latitude, and present 
the results in Sect. 15.31 

5.2. The 2D calculation of the emerging spectrum 

Ideally, solving a full set of 2D radiative transfer and statisti- 
cal equilibrium equations would be required to analyze a non- 



5.3. Model fits, results, and constraints on the wind 
asymmetry 

The 2D model spectra were computed following the physical pa- 
rameters of the best CMFGEN spherical model shown in Tables 
[2] and [3] We varied the density enhancement and velocity field 
as a function of latitude, and Table [5] summarizes the properties 
of the 2D models. Figure [9] shows the latitudinal changes in the 
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Table 5. Summary of the paramet ers used to compute the spec- 
trum in 2D geometry using the Busche & Hilliej (120051) code. 
Note that the density contrast between equator and pole is given 
by 1 + a for oblate models, and 1/(1 + a) for prolate models. 



Model 


Density 
profile 


a 


b 


Velocity 
changes? 


VB\ 


VB2 


VB3 


1 


oblate 


1 


2 


no 








2 


prolate 


1 


2 


no 








3 


oblate 


3 


2 


no 








4 


oblate 


3 


5 


no 








5 


oblate 


7 


2 


no 








6 


oblate 


7 


5 


no 








7 


oblate 


15 


2 


no 








8 


none 






yes 


0.2 


0.8 


4 


9 


oblate 


3 


2 


yes 


0.8 


2.2 


4 



density enhancement (normalized to the best spherical model) 
for the different models. We analyzed the effects of latitudinal 
changes in the density by comparing the 2D model spectra with 
the observed line profile of Hei A 5876, which is very sensitive 
to latitudinal changes, and has the most deviating line profile in 
the spherical model. This line is also isolated, minimizing the 
errors due to blending. Similar effects are seen in other low ion- 
ization lines of He i, C m, and N m. 

Figure \W\ displays the observed spectrum of the qWR star 
around 5800 A compared with different 2D model spectra, la- 
beled as in Table [5] All models were computed for a viewing 
angle of i = 0.77° (Paper I), and the best spherical CMFGEN 
model is overplotted on each panel for comparison. 

The immediate conclusion from Fig. [10] is that the wind of 
the qWR star cannot be prolate (model 2), as the fit to the He i 
line becomes even worse than the fit from the spherical model. 
This happens because, as the system is viewed pole-on, a prolate 
wind implies an enhanced density in the polar regions compared 
to the spherical model, producing stronger P-Cygni absorption 
and stronger emission at high velocities - exactly what is seen in 
the spectrum of model 2. 

We can also rule out that the narrow Hei line profile is 
only due to changes in the wind terminal velocity as a func- 
tion of latitude, as is illustrated by model 8. This model assumes 
a slow equatorial wind (voo=70 kms~') and a fast polar wind 
(voo=700 kms -1 ). Although the fit to the emission component 
is improved, the P-Cygni absorption is still very strong, since 
the density enhancement has not been changed as a function of 
latitude. 

Therefore, it seems clear that a density enhancement in the 
equatorial region is required to reproduce the observed emission, 
and a density depletion along the polar regions is required to 
reproduce the lack of a P-Cygni absorption in Hei A 5876. The 
minimum density contrast equatonpole that fits the observations 
is 8:1, but we cannot rule out larger enhancements, such as 16:1. 
However, we have to keep in mind that larger enhancements will 
likely change the ionization structure of the wind as a function of 
latitude, which is not accounted in the current calculations. The 
parameter b, which describes how fast the density changes from 
equator to pole, is quite difficult to be constrained. A model with 
a density contrast of 8:1 and b — 5 (model 6) produces almost 
the same He I line profile as a model with a density contrast 16:1 
and b — 2 (model 7). 

We also examined changes due to both an oblate density en- 
hancement and changes in the wind terminal velocity as a func- 
tion of latitude. However, there is no clear evidence from the 
fits to the optical lines that a higher velocity polar wind is re- 
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Fig. 10. Comp arison between the line profiles computed with 
the 2D code of iBusche & Hillieri (120051) (blue line) with the ob- 
servations of HD 45166 around Hei A 5876 (black line). For 
each panel the models are labelled as in Table |5j followed by 
the density contrast equator:pole, and the parameter b. The best 
spherical-symmetric CMFGEN model (red dashed line) is re- 
peated in each panel for comparison (see text for further details). 
The weak dip in the observed Hei A 5876 line profile is likely 
due to incomplete subtraction of the secondary star spectrum. 



quired, although it cannot be ruled out. If there is such change in 
the wind terminal velocity, the minimum density contrast from 
equator to pole is reduced to 4: 1 (see model 9 in Fig.fTOt. Clearly, 
the largest effect of a fast polar wind would be seen in the ultra- 
violet resonance lines. However, we refrained from analyzing 
the archival IUE data simultaneously to our optical dataset due 
to the high variability of HD 45166 and the long time lag of 
15-20 years between the IUE ultraviolet observations and our 
optical data. We acknowledge the wealth of information and the 
fundamental importance of the UV observations to constrain the 
nature of HD 45 166, and a forthcoming paper will be devoted to 
a detailed analysis of the UV spectrum. 

It is worthwhile to stress that there is clear observational ev- 
idence of a high velocity outflow in HD 45166 through the anal- 
ysis of the discrete abso rption components (DAC) in the reso- 
nance lines of HD 45166 dWillis et alJ ll989). However, it is also 
clear that marked variability is present, and at some epochs the 
UV resonance lines (e.g. Civ AA 1548-1550) do not show any 
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Fig. 9. Density contrast of the HD 45 1 66 models with latitude-dependent density distributions (normalized to the density distribution 
of the spherical model), as a function of latitude. Each model is labeled as in Table |5J model 9 has the same density distribution 
as model 4, and is not shown. The equator of the star is at [x, 0], while the pole is located at [0, v]. The colorbar in the right has a 
linear color scale, which we adopted the same for all models, to highlight the different density variations as a function of latitude. 
The maximum value in the colorbar is 1.8, which corresponds to the maximum value found in model 5. 



high velocity absorption component dWillis et"al] Il989t) . Two 
questions need to be answered. 

1. Does the velocity of the DACs reflect the true wind termi- 
nal velocity along the polar direction? While in OB stars the 
maximum velocity of the DACs tend to mim ic the wind ter- 
minal velocity (e.g. lPrinia & Howarfhl ll986). further work is 
required to establish whether the same holds for HD 45166. 

2. Do the DACs reflect an intrinsic property of the qWR and its 
steady wind or is the presence of a close companion some- 
how related to the variability seen in the UV lines? 

Contemporaneous ultraviolet and optical spectroscopy are re- 
quired to answer these questions, and we will come back to them 
in Paper III. 

We have not considered in this work latitudinal changes 
in the wind ionization structure derived from the spherical 
CMFGEN modeling (Sect. 03]), which likely occur in the case of 
large density changes. The inclusion of this effect would likely 
reduce the density contrast determined from our analysis, but the 
main conclusions reached in this work would remain valid. This 
is because a density contrast is needed to change the ionization 
structure of the wind, since the radiation field of the B7 V com- 
panion is too weak. The exact treatment of the effects of the lati- 
tudinal changes in the ionization structure of the wind is beyond 
the scope of this paper. 

6. The Temporal Variance Spectrum 

We performed the Temporal Variance Spectrum (TVS) analysis 
in order to study the characteristics of the emission line profiles. 
In this procedure, the temporal variance is calculated, for each 
wavelength pixel, from the residuals of the continuum normal- 
ized spectra and the average spect rum. For further deta ils and 
discussion on the TVS method, see lFullerton et alj d!996l) . 

We calculate the TVS for the strongest spectral lines, such 
as Hen A 4686, Hei A 5876, and Civ AA 5801-5811 (Fig.QTJ- 
The TVS contain a large amount of information that, in spite of 
not being always of immediate interpretation, certainly can be 
helpful in understanding the wind structure. Initially, one should 
notice that, as well as in the line profiles, also in the TVS the 
spectrum is quite different when we compare, for instance, He n 
A 4686 and C iv AA 5801-5812. At the wings of He n A 4686, the 
TVS is close to a Lorentzian profile. At smaller velocities, how- 
ever, the TVS profile of the line deviates more and more from 
a Lorentzian profile until it reaches a local minimum. Assuming 



that the Lorentzian profile is a signature of optically-thin line 
emission, the wind seems to become transparent for He n line 
photons at velocities larger than 270 kms" 1 . The correspond- 
ing radius is r — 10/?*, according to the velocity law determined 
using the best CMFGEN model. The CMFGEN model itself pre- 
dicts that r(He n4686) < 1 for r > 15 fl*(see Fig.!). 

On the other hand, the TVS of the Hei A 5876 line presents 
a peculiar behavior. While the intensity profile is asymmetric 
to the red, the TVS profile is asymmetric to the blue. This dis- 
crepancy is quite strong and its interpretation is not trivial. The 
stronger H lines of He and H/3 are blended to the He n lines and, 
therefore, their profile and interpretation is more complex. 

A noticeable feature in the TVS profiles is the presence of 
a central dip in the Hen and Hei lines. In the case of Hen, 
this dip is centered at a velocity of +35.3 + 0.6 kms . We do 
not have an interpretation for this velocity. In the TVS profile 
of Hei A 5876, however, the dip is at a different velocity of 
+4.1 ± 0.1 kms -1 , which is consistent with the systemic veloc- 
ity of +4.5 ± 0.2 kms 1 (Paper I). On the other hand, in the 
intensity spectrum of He i A 5876, the absorption is centered at 
+9.2 ± 0.3 km s . Interestingly, this is remarkably similar to the 
photospheric velocity found for the absorption lines of the sec- 
ondary star. Therefore, we suggest that such dip in the TVS of 
He i A 5876 is due to the orbital motion of the secondary star. 

Most of the intensity line profiles present in the spectrum of 
HD 45166 are quite symmetric. The lines of Civ A A 5801-581 1, 
for instance, are very well fitted by symmetrical Lorentzian pro- 
files, as already mentioned in Paper I. However, in the TVS the 
asymmetries are much more evident. For the lines of Nm AA 
4634^1640, C in AA 4647^1650 and Crv A 4658, as well as for 
C iv AA 5801-5812 (see Fig. QTJ, the TVS profiles are very dif- 
ferent from the intensity profiles of the respective lines as well 
as from the TVS of the He n A 4686 line. While in the TVS of the 
latter the blue wing is 30% more intense than the red wing, in the 
lines of C and N the opposite happens. For instance, in the TVS 
of C iv AA 5801-58 12, the red wings are 3 times stronger than the 
blue wings (see Fig. ITTb . We interpret this phenomenon as evi- 
dence for temporal changes in the gas density of the equatorial 
wind region. This produces a variable optical depth (due to ab- 
sorption or electron scattering) and primarily affects the red part 
of the profile, which is formed in the receding side of the wind. 
The nature of the variability in the gas density is not clear, and 
a detailed interpretation is beyond the scope of this paper. Two 
possible interpretations are variability in the (equatorial) mass- 
loss rate, or debris from a hypothetical, variable mass transfer 
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Table 6. cr/I for the strongest lines present in HD 45166. Both 
<x and / were measured at the centroid of the lines. 
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Fig. 11. Average intensity spectrum and TVS of HD 45166 for 
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and dot-dashed horizontal lines, respectively. The vertical dotted 
lines mark the position of the zero velocity for the relevant line 
presented in each panel. See text for details. 



from the secondary star. The latter hypothesis will be examined 
in Paper III. 

A characteristic indicator of each line is the ratio between its 
variance and the intensity, cr/I. Table|6]presents this ratio for the 
strongest lines of HD 45 166. Q It is interesting to notice that Hei 
lines have the highest cr/I values (2.6%), followed by Hi+Hen 
lines, which are estimated to be about 2.5%. The higher variabil- 
ity of low ionization lines can be interprete d in the context of a 
Wind Compression Zone (WCZ) scenario (Ign ace et al.| [l996). 
in which a lower ionization region is produced at the equator of 
the system (see Sect. 17.2b . Variability in the WCZ may cause a 
higher cr/I ratio for the He i and H i lines when compared to the 
higher ionization lines of He n or C iv. It is not clear whether the 
same scenario can explain the variability detected in the p hoto- 
spheric UV lines of Fe v, He II, and N iv dWillis et al.ll 19891) . 



7. Discussion 

7.1. Gravitational redshift and the line formation regions 

The gravitational redshift is measurable in the CNO emission 
lines of HD 45166 (see Paper I) . The CNO III lines have an av- 
erage radial velocity of 5.0 + 0.8 kms 1 (Tabled. This is com- 
patible with the systemic velocity y = 4.5 + 0.2 km s _1 (Paper I). 
However, the N iv-N v lines have an average radial velocity of 
8.3 + 1.2 kms 1 (Tabled These higher-ionization lines should 
be emitted in the wind inner layers and, therefore, should be red- 
shifted when compared to the lines of lower ionization. Indeed, 



1 The telluric lines in absorption have a much higher TVS in our data 
(28%), and, therefore, those lines can be easily identified and differ- 
entiated from the stellar lines. On the other hand, the interstellar lines 
should not appear in the spectrum of the TVS. The TVS can be a useful 
method to distinguish between the several types of lines present in a 
rich emission and absorption line spectrum such as in HD 45 166. 



Ion 


Line 


cr/I (%) 


Hi+Hen 


4861 


2.0 


Hi+Hen 


6563 


2.2 


He i 


5876 


2.6 


He i 


6678 


2.1 


Hen 


4541 


1.6 


Hen 


4686 


1.5 


Hen 


6527 


1.4 


Hen 


6683 


1.3 


Civ 


5801 


1.6 


Civ 


5812 


1.5 


Cm 


4647 


1.9 


Nm 


4510 


1.1 


Nm 


4640 


1.9 



Table 7. Average FHWM and radial velocity of the spectral lines 
for the most abundant ions present in HD 45166. 



Ion 


< FWHM > 


< V r > 




(km s" 1 ) 


(kms" 1 ) 


Civ 


142 ± 19 


4.4 ± 0.4 


Cm 


91 ±6 


5.4 ± 1.5 


Nv 


79 ± 20 


7.3 ± 1.6 


Niv 


82+11 


8.6+ 1.5 


Nm 


96 ±4 


4.0 + 1.1 


Om 


113 ± 9 


5.0+ 1.6 


< CNO in > 


102 ± 4 


5.0 ± 0.8 


< N v + N IV > 


78 + 10 


8.3 ± 1.2 



the observed velocity shift of those more ionized species, com- 
pared to the less ionized ones, is 3.3 + 1.2 kms -1 . This means 
that, taking the systemic velocity as a reference, the CNO III 
lines should be emitted at distances larger than ~ 10 as- 
suming the velocity law derived in Section 14.31 and M q wR = 
4.2 M (Paper I). Interestingly, those distances are comparable 
to the size of the Roche lobe of the qWR star. In contrast, the 
gravitational redshift of the Niv-Nv lines indicates that they 
should be emitted close to the hot star, at distances of the order 
of 1.5 -2.0/?*. 

The above values are in reasonable agreement with the line 
formation regions presented in Fig. [8] taking into account that 
the lines actually form in an extended wind region. One should 
also consider that the emission lines can be a result of pure re- 
combination (e. g. Civ A 4658) or, at least partially, due to con- 
tinuum fluorescence (e.g. Civ AA 5801-5812). This is shown 
in Fig.Q/JJ where the line formation region of Civ AA 5801 is 
compared to the radial variation of the density of C 3+ and C 4+ 
ions. It can be seen that roughly 50% of the line emission come 
from an inner region where the dominant ionization stage of car- 
bon is C 4+ , suggesting that the line is formed by recombination. 
However, a significant fraction of the line formation occurs un- 
der very low C 4+ density conditions, making it very unlikely that 
recombination is the main mechanism of line emission in those 
regions. Instead, we suggest that the line is emitted due to con- 
tinuum fluorescence of the C 3+ ions. 

This discussion illustrates that the physical conditions found 
in HD 45166 is far too complex to be properly described by a 
spherical wind model. Therefore, understanding the line forma- 
tion processes and determining the line formation regions pro- 
vide invaluable information on the nature of HD 45 166. 
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log (r/R.) 



Fig. 12. Upper panel: radial density profile of the C 3+ and C 4+ 
ions along the wind of HD 45166. Lower panel: line formation 
region of Civ AA 5801. The quantity £ is related to the EW of 
the line (following |miiieJll989t) as EW = j£° £(r)d(log r) . 



7.2. Rotation, wind asymmetry, and formation of a Wind 
Compression Zone 

Depending on the rotation velocity of the star and on the 
velocity law of the wind, a wind compression disc (WCD, 
Biorkman & Cassinelli 1993), or a wind compression zone 
(WCZ. Hgnace et all l996). can be formed. The most fundamen- 
tal parameter to consider is the ratio of the rotational velocity 
of the star (v 10t ) to the critical break-up velocity (v cl -i t ). Since 
HD 45166 is seen almost pole-on (i = 0.77°, Paper I), it is 
impossible to detect the rotational broadening of lines formed 
close to the photos phere, and to fur t her de termine the rotational 
velocity using the lBusche & Hillien (120051) code, similar to what 
has been done, for instanc e, for the luminous blue variable AG 
Carinae dGroh et al.ll2006l) . 

Nevertheless, the effects due to the rotation of HD 45 166 can 
be quantified in terms of the ratio of the critical period for break- 
up (P cr i t ) to the rotational period (P rot ), 

Vrot/Vcrit = Pcnt/Ptot , (4) 

where Vcrit = (GM+/R+) 112 . For values of v ro t/Vcrit clo se to 1 
a WCD can be formed (Biorkman & Cassinellil 11993b . while 
for moderate rota tional velocities a WCZ can be produced 
dlgnace et al.ll 1996b . 

Using the stellar parameters determined for HD 45166 (see 
Sect. HJ, we have P cl -i t ^ 1 .0 h. Therefore, unless the rotational 
period is close to 1 hour, the rotation is too slow to cause the 
formation of a WCD. However, the formation of a WCZ may 
still be possible, and might explain the degree of asymmetry de- 
tected in the wind of HD 45166 (see Sect. [5}. In Paper III we will 
show evidence of a period of 2.4 hours, which might be related 
to rotational modulation. 

The difference between the two scenarios outlined above is 
that, in the case of a WCD, the formation of a stationary shock 
heats up the disc, and compresses i t to a thin sheet with high den- 
sity dBjorkman & Cassinellilll993l) . In the case of a WCZ, this 
shock does not exist and the comp ression area is much thicker 
and less dense (Ignace et al.lll996b . While for a WCD the ratio 
between equatorial and polar density is of the order of ~ 1000, 



for a WCZ with the parameters of HD 45 1 66 this ratio is ~ 5 — 10. 
The value of the density contrast derived from our 2D modeling 
is 4-8 (see Sect. |5j, which is of the same order as predicted by a 
WCZ scenario. 

In addition, in the case of a WCD, the shock formed in the 
equatorial region should be hot enough to be detected in soft X- 
rays. However, HD 45166 was not detect ed by the ROSAT X-ray 
satell ite. Superionization due to X-rays dBjorkman & Cassinellil 
1 1993b should also produce enhanced N v and O vi emission near 
the equator, which is not observed. In contrast, on the equator 
there is enhanced emission of low ionization species such as He i 
lines (see Sect. [5]l, which argues against the presence of a WCD. 

However, there are at least two caveats which might change 
the WCZ/WCD scenario and have to be taken into account 
in future works of HD 45166. First, non -radial line forces 
might play an important role in shaping the 2D wind structure 
(iGavlev & Owo cki 2000). Secondly, if the star is rapidly rotat- 
ing, the effe ctive temperature is also going to vary as a function 
of latitude (Ivon Zeipellli924h . and the 2D wind structure will 
likely be affected. 

8. Conclusions 

In the following, we present the main conclusions of this paper. 

1. We used the radiative transfer code CMFGEN to quantita- 
tively analyze the wind of the qWR star HD 45166, and to 
constrain its fundamental parameters. Comparing the model 
spectrum with high-resolution optical observations, it was 
possible to determine the effective temperature, radius, lu- 
minosity, mass-loss rate, wind terminal velocity, and sur- 
face abundances. The temperature of the hot star is T e g = 
50000 + 2000 K, the radius is R phot = 1 .00 R , and the lumi- 
nosity is log(L/L ) = 3.75 + 0.08. 

2. The wind parameters are quite anomalous. The mass loss 
rate determined is M = 2.2 x 10~ 7 M G yr _I , which is about 
4 time s higher than the value proposed by Ivan Blerkoml 
(1978). Analyzing the optical lines, we determined the ter- 
minal velocity of the wind t o be v ro =42 5 kms " 1 , and not 
1200 kms~'as measured by Will is et aD dl989h . since the 
optical spectrum shows no evidence of such a high-velocity 
wind. While for O-type and WR stars Voo/v esc > 1-5, in the 
case of HD 45166 it is much only 0.32. In addition, the effi- 
ciency of momentum transfer from the radiation field to the 
wind is r\ — 0.74, which is at least a factor of 4 lower than 
what is found for WR stars. 

3. The star is helium-rich, with Nn/Nne — 2.0, obtained by fit- 
ting simultaneously the He n lines which are unblended, such 
as He ii A 4686, and those blended with hydrogen Balmer 
lines, such as He n+Ha A 6560. The He content is consistent 
with the analysis of the He II Pickering lines, which shows 
that a significant abundance of hydrogen exists in the wind. 
The helium content of HD 45166 is actually 4 times h igher 
than the abundance determined bv lvan Blerkoml dl978 ). The 
CNO abundances are quite peculiar, and very different from 
those of the Sun, of central stars of planetary nebulae, as well 
as from WN stars. 

4. The wind of the qWR star is not spherically symmetric. 
Analyzing the observed line profile by comparison with 
model spectra computed in 2D geometry, we were able to 
infer the presence of an oblate wind, with a minimum den- 
sity contrast of 8: 1 from equator to pole. Taking into account 
our assumptions, higher density contrasts, for instance 16:1, 
cannot be discarded since they yield similar line profiles. 
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5. We also analyzed whether the wind terminal velocity 
changes as a function of latitude, but there is no clear evi- 
dence that this is the case. If we include a fast polar wind 
(v M =1300 kms -1 ) and a slower equatorial wind (v oo =300 
km s _I ), the observed line profiles can be fitted with a smaller 
density contrast from equator to pole of 4: 1 . Therefore, such 
possibility might be a promising way to reconcile the values 
of v M obtained from optical and UV analyses. 

6. We applied the Temporal Variance Spectrum (TVS) tech- 
nique to the strongest emission lines of HD 45166. The 
TVS profiles of the lines is quite complex. The ratio 
sigma/intensity varies from 1.5% for Hen lines up to 2.5% 
for those of Hei. Both Hei and Hen lines also present a 
central dip in their TVS, in spite of this signature not being 
present in the intensity spectrum of He n. 

7. The TVS of some emission lines, especially C iv, show that 
variable absorption or scattering occurs near the equator; as 
a consequence, the red wings are more variable than the blue 
wings of the same line. The cause of this variability could 
be blobs of gas moving out in the wind compression zone or 
debris from variable mass transfer from the secondary star. 

8. Differential gravitational redshift seems to be noticeable in 
the velocities of high-ionization lines when compared to the 
velocities of low-ionization lines. 

9. The latitudinal changes in the density and velocity might be 
explained by the f ormation of a Wind Compression Zone 
(Igna ceet al.| [l996). which would imply that the qWR star 
has a relatively high rotational velocity. 
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Appendix A: The H/He abundance and the He II 
Pickering decrement 



The hydrogen abundance determined by Ivan Blerkoml (Il978h 
is Nn/Nne ~ 5.3 (by number), which is about 4 times higher 
than the value derived through the detailed radiative transfer 
modeling using CMFGEN (N H /N He =2.0 by number, see Sect. 
4.41 ). This difference is significant, and we will analyze the Hen 
Pickering decrement to verify this discrepancy. 

Although the hot star of HD 45 166 is He-rich ( 65% in mass, 
at least in the wind), the presence of hydrogen is an important 
as pect to be carefully consid ered in the analysis. As mentioned 
bv lWillis& Stickland! (119831) . jud ging from the He n Pickering 
decrement, hydrogen is definitely present. To check this, we 
compared our observations with the expected theoretical values. 
For practical reasons we normalized the predicted and the ob- 
served values to 1 for Hen A 5411. Figure [ATI shows the com- 
parison between the two series, and we clearly see the typical 
oscillation pattern that appears when hydrogen is present. It is 
also possible to infer that the intensities of the lines that coin- 
cide with hydrogen Balmer lines are about twice as strong as the 
expected ones if there was no Balmer emission. One way of ana- 
lyzing the presence of hydrogen in a helium-rich wind spectrum, 
following lOliveira & Steiner! d2004l) . is by defining the Pickering 
parameter, p: 



7(4859 + 4861) 
[7(4541)./(5411)]° 



(A.l) 



Using the predicted theoretical intensities that one expects 
for a pure He n spectrum, calcul ated for a low-density gas, 
p = 0.97 ± 0.01 (IOsterbrocklll989l) . Any value of p larger than 
1 would mean hydrogen is present in the wind. In the case of 
HD 45166, the measured value is p 2.95. In case both H and 
He lines were optically thin, it would be easy to determine the 
relative abundance between the two species. In this case, given 
the sp ecific emissivity of the respective transitions dOsterbrockl 
1989), 



= p-l 



(A.2) 



N(H+) 

N(He ++ ) 

ForHD45166,onegets/V(// + )/A?(i/e ++ ) 1.95. The largest 
uncertainty in this value is the hypothesis that all of the involved 
lines are optically thin. This may be a good approximation at the 
wings of the li nes but not at th eir cores. For an optically thick 
case we have (IContietal.lll983h : 



N(H+) 



1 



N(He ++ ) ' 
and thus we obtain N(H + )/N(He ++ ) 



(A3) 



4.06. 
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Upper quantum number 

Fig.A.l. Intensity of the Hen lines of the Pickering series as a 
function of the upper quantum number, normalized such as that 
I(He ii A 54 1 1 ) = 1 . The open squares represent the observed val- 
ues of HD 45166, while the filled circles cor respond to the the - 
oretical prediction for an op tically-fhin gas (lQsterbrock| [l989). 
using tabulated values from lOliveira & S teiner (120041) . 



Appendix B: The CNO emission line table for 
HD 45166 
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Table B.l. Equivalent width of the CNO emission lines in HD 45166 obtained from the observations. References: Willis & Stickland (1983) for 
ultraviolet lines (1200-3300 A), and this work for optical lines (3700-9000 A). Key: y=line is present; y?=doubtful identification due to low S/N; 
?=doubtful identification due to presence of a telluric line; bl=blended with other lines; n=line is definitely not present. 
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